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Abstract We present a grid of stellar models, obtained with 
the CLES evolutionary code, following the specification of 
ESTA-Taskl, and the corresponfing seismic properties, com- 
puted with the LOSC code. We provide a complete descrip- 
tion of the corresponding files that will be available on the 
ESTA web-pages. 
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1 Introduction 

The preparation of the CoRoT mission brought up the need 
of a reference grid of stellar models in order to locate the 
possible CoRoT targets in the Hertzsprung-Russell (HR) dia- 
gram, and to allow a first interpretation of the forthcoming 
CoRoT-data. For this purpose, several grids of stellar models 
were computed according to the ESTA-Taskl specifications. 
In this paper we present the grids of stellar models we com- 
puted with the code CLES (IScuflaire et all l2007ah . as well 
as the seismic propert ies we derived by usin g the adiabatic 
oscillation code LOSC (IScuflaire et al.Ll2007bl) for models on 
the main sequence (MS) or close to it. 

In sections |2] and [3] we describe, respectively, the input 
physics and the set of stellar parameters we used. The evolu- 
tionary tracks, stellar models and oscillation frequency data 
can be found on the ESTA Web site. The description of all 
this material is given in Sect|4]for the evolution and struc- 
ture models, and in Sect|5]for the seismic properties of the 
stellar models. Finally, in Sect.|6]we summarize the type and 
amount of available files. 
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2 Input physics 

For three out of four stellar model grids we adopted the ref- 
erence input physics, as well as the physical and astronomi- 
cal const ants specified fo r ESTAT ask 1 comparisons (see also 
Sect.5 in lLebreton et al.l [200731) . We stress that these spec- 
ifications do not always coincide with those of the standard 
version of CLES (see Scuflaire et al., 2007a). 

- Equation of State (EoS) - O PAL 2001 equation of State 
(iRogers & Navfonovl l2002h available in the OPAL Web 
site. The tabulated values of Cy have been replaced with 
the values derived from the tabulated values of P, Fi , Xt, 
Xp- For each p, T, X and Z we interpolate in the tables 
by using our own interpolation routi ne. 

- Opacities - OPAL96 ( Iglesias & Ro gersLll99d) . The opac- 
ity tables used in CLES have been calc ulated on-line for 
the standard GN93 tGrevesse & Noelslll993h metal mix- 
ture, and using the smoothing routine available in the 
OPAL Web site to build the final opacity tables. The l ow 
temperature opacity tables bvlAlex ander & FergusonI (Il994) 



are also smoothly added (see Scuflai re et al 



2007ah .and 



conductive opacities were not included in these compu- 
tations. As for EoS, we interpolate in p, T, X and Z by 
using our own interpolation routines. 

- Nuclear network - We choose the same nuclear net- 
work as in model comparison, that is, basic pp chain and 
CNO cycle reactions up to the ^^0(p, a) '^N. We adopted 
the nuclear reaction rates from the analytical formulae 
provided by the NACRE compilation dAngulo et alll 19991) . 
including that one for the '^N(p, y)'^0 , and the weak 
screening factors from ISalpetej (119541) . The electronic 
density, taking part in the screening factor computation, 
was estimated assuming full ionization of chemical ele- 
ments. 

- Chemical compositio n - We adopted the meta l distri- 
bution provided by the lGrevesse & Noelsl (Il993h (there- 
after GN93) solar mixture. Only the abundances of light 
elements (Li, Be and B) were modified to be consistent 
wi th those used in CESAM. T he isotopic ratios are those 
of lAnders & Grevessel (Il989h except for the ^H/'H and 
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Fig. 1 Theoretical HR diagram for the evolution phase comprised between the initial homogeneous model and that with a hydrogen mass fraction 
at the center reduced by 1% with respect to the initial model (that is, Xc = 0.693). Evolutionary tracks are labeled by the corresponding stellar 
mass. Solid lines correspond to grey models and dashed ones to models with Kurucz boundary conditions. 



^He/^He ratios for which we took those of lGautier & Morell - 
([1993). 

Convection - All the evolutionary tracks were computed 
wi th the classic al mixing length treatment of c onvection 

Bohm-Vitense C1958.) and the formulation of lHenvev et al.l 
(Il965h for optically thin regions. The mixing length pa- 
rameter aMLT was fixed at 1 .6 for all the grids. Note that 
the latter is the value chosen for ESTA comparisons, but 
it does not correspond to the value derived from a solar 
calibration. 

Overshooting. All the models, regardless of their mass, 
were computed without convective core overshooting. 



Atmosphere. ESTA-Taskl specification requires the sur- 
face boundary conditions to be provided by integration 
of a grey atmosphere following the Eddington's ^(t) 
law (grey models, thereafter). Three of the four stellar 
model grids were computed following the ESTA speci- 
fication, while for the fourth one, temperature and den- 
sity at r = Tf-ff were obtained from Kurucz's atmosphere 
models (.Kurucz„,199&) . 
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Fig. 2 As Fig. [T] for the evolutionary phase of main sequence and post-main sequence. Open circles indicate the location of models that has 
burned, at the center, half of the initial hydrogen. Solid circles correspond to models with only 0.1% of hydrogen at the center. 



3 Stellar parameters 



We provide evolutionary tracks for masses from 0.8 to 8 Mq, 
with a mass step AM=0.l, from 0.8 to 1.6 Mq, AM=0.2 up 
to 4.0 Mq, and ZiM=0.5 from 4.0 to 8.0 Mq. 

The initial hydrogen mass fraction is X = 0.70 in all 
grids. Three different values of the metal mass fraction Z 
(0.02, 0.01 and 0.006) are available for grids of grey models, 
while only Z = 0.02 has been considered in the grid com- 
puted with Kurucz atmosphere boundary conditions. 



4 Data 

4.1 Evolutionary tracks 

For each stellar parameter (M,X,Z) we follow the PMS evo- 
lution from the Hayashi track, and the calculation ends either 
when the temperature at the stellar center is high enough 
to burn He, or when the age of the model is larger than 
20 Gyr. We recall that the computation of these models does 
not include conductive opacities and that these ones should 
be taken into account in modeling low mass stars up to the 
helium ignition. The grids of these evolutionary tracks are 
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available from the ESTA Web sit43 and consist of two HRD- 
files, whose name contains the values of the stellar parame- 
ters for which they were computed. For instance, "ml .00Z0.02X0.70- 
HRDl .txt" and "ml .00Z0.02X0.70-HRD2.txt" contain the evo- 
lutionary tracks for a 1 .00 M© star with initial chemical com- 
position X = 0.70 and Z = 0.02. The first four lines in HRDl 
and HRD2-files provide information about the stellar param- 
eters and the input physics used to build the model. Since 
most of them have already been described in previous sec- 
tion, we only mention that fg and ft indicate the factor by 
which the mesh and time step were multiplied. Both, fg and 
ft were set equal to 0.5, that is, the number of mesh points 
and the temporal steps were dou bled with respec t to the stan- 
dard value used in CLES ( S cuflaire et al.l l2007bh . The quan- 
tities for which the temporal evolution were tabulated in the 
HRDs files are the following: 
HRDl files: 

- Column 1 - Stellar mass in units of Mq. 

- Column 2 - Stellar luminosity (log(L/L0)). 

- Column 3 - Decimal logarithm of the effective temper- 
ature in K (log Teff). 

- Column 4 - Stellar radius in units of solar radius {R /Rq). 

- Column 5 - Stellar age in Myr. 

- Column 6 - Hydrogen mass fraction at the center. 

- Column 7 - logg, with g, the surface gravity . 

- Column 8 - Number of convective -radiative boundaries. 

- Column 9 - Codification of the of boundary type: "2" 
for boundary from convective to radiative region ;and 
"1" for boundaries from radiative to convective one, go- 
ing from the stellar center towards the surface. 

- Column 10-15 - Relative radius of the six first convective- 
radiative boundaries (from the center outwards). 

HRD2 files: 

- Column 1 - Number of model. 

- Column 2 - Stellar age in Myr. 

- Column 3 - Stellar radius in solar radius (R/Rq). 

- Column 4 - Stellar luminosity in units of solar luminos- 
ity (L/Lo). 

- Column 5 - Effective temperature in K. 

- Column 6 - Temperature at the center in units of lO' K. 

- Column 7 - Central density pc in g.cm^^. 

- Column 8 - Hydrogen mass fraction at the center (X^.). 

- Column 9 - Central helium mass fraction (Yc). 

- Column 10 - Number of convective regions. 

- Column 11 - Relative mass of the convective core (mcc/Af* 

- Column 12 - Relative radius of the base of the outer 
convective region (7?ce/^*)- 

- Column 13 - Mass of the helium core in Mp,. The he- 
lium core is defined by the mass of the region with a 
hydrogen mass fraction smaller than 10~^. 

The logTeff-logL/L© diagram for all the models with a 
metal mass fraction Z = 0.02 have been plotted in Fig. [T] 
for the PMS evolution, and in Fig. |2] for the MS and Post-MS 
evolutionary phases. For each mass there is a pair of curves, 

' http://www.astro.up.pt/corot/models 
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Fig. 3 Relative mass of the convective core as a function of the stellar 
mass for models that have burned the half of their initial H content at 
the center. 




(MJ 

Fig. 4 Relative radius of the outer envelope bottom as a function of the 
WUar mass for the same models as in Fig. [3] 



one was computed with Eddington's law as boundary con- 
ditions at the photosphere (solid lines) and the other used 
Kurucz's atmosphere models (dashed lines). It is interesting 
to note that, as it concerns the theoretical HR diagram, the 
atmosphere type affects only models with effective tempera- 
ture lower than 6300 K. The effects on the seismic properties 
will be discussed in Sect. [5] 

In Fig. |2] we marked the model locations at the middle 
of their MS (central hydrogen mass fraction X^ = 0.35) by 
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Fig. 5 Relative mass of the helium core at the end of the MS phase as 
a function of the stellar mass. 



open circles and that of models with a central hydrogen mass 
fraction ~ 0.001 by filled ones. Defining the beginning 
of MS by the model with a central mass fraction of hydro- 
gen decreased by 1 % with respect to the initial one (that is 
Xc = 0.693 for the present computations), and its end by that 
withXc ~ 0.001, the MS-lifetime varies from 5 x 10^ yrs for 
8 Mq to more than 13 Gyr for 0.9 Mq. The number of com- 
puted stellar models that span the main sequence phase of 
evolution is of the order of 160 for each (M,Z,Z). 

The dependence on the stellar mass of the convective 
core mass, for models in the middle of MS is shown in Fig. [51 
and that of the convective envelope bottom in Fig.|4l Finally, 
in Fig.|5]we show the mass of the He core at the end of the 
MS as a function of the stellar mass. 



4.2 Stellar models 

In addition to the evolutionary tracks, a selection of internal 
structure models for each {M,X,Z) are also available : 

- ZAMS model - defined as the model for which the ratio 
of gravitational to nuclear luminosity is lower than 2%. 

- IMS model - as the model which has burned 1 % of H 
(Xc =i 0.693). 

- MS model - as the model with half of the initial H at the 
center (X^ =i 0.35) 

- TAMS model - as the model with X^ ~ 0.01 

- One every five models spanning the MS phase of evolu- 
tion: from Xc = 0.693 to Xc = 0.001. That means about 
30 models for each set of stellar parameters. 



ml.00Z0.02X0.70-####.gong, where "####" can be "ZAMS", 
"IMS", "MS", "TAMS", or a number corresponding to the 
number of the model. The format adopted for these files is 
the FGONG0 one with the first twenty-five variables. The 
number of mesh points through the stellar interior is of the 
order of 2200, and 100 additional mesh points are used to de- 
scribe the stellar atmosphere from an optical depth T = 2/3 
outwards, up to T = 10^^. 

Among the physical quantities stored in FGONG files 
there are the mass fraction of the most important contrib- 
utors to nuclear energy generation. For illustration, we show 
in Figs.|6]the chemical profile of key elements such as ^He, 
'^C, and '^N inside the stellar models in the mid-MS. 



5 Seismic properties 

For each equilibrium model in FGONG format, we also pro- 
vide a file containing the properties of radial and non ra- 
dial oscillation modes with spherical degree £ = I, 2, and 
3. These properties wer e computed by using the code LOSC 
dScuflaire at all [2007b) with the standard surface boundary 
condition (regularity of solution when P = at the surface, 
that is, 5P/P+ (4 + (0^)5r/r = 0). The angular frequencies 
of the computed oscillation modes cover the domain given 
by 0.3-50 times the dynamical time (Tdym = y^/GM). 

The file name is the same as that of the equilibrium model 
but with the extension ".freq"; the content is the following: 

- Column 1 - spherical angular degree I. 

- Column 2 - radial order n. 

- Column 3 - dimensionless frequency ft) = a • Td. 

- Column 4 - angular frequency a 

- Column 5 - frequency V = a/2 tt in Hz. 

- Column 6 - period P in seconds. 

- Column 7 - /3„^ = 1 — Cni gives the frequency shift gen- 
erated bY_soIid rotation, where C„^ is the Ledoux's con- 
stant (L edouxL[l95lh . 

- Column 8 - fraction of the kinetic energy associated to 
the radial component of the motion ev — E^^^y /£^kin 

- Column 9 - = (x) is the mean value of x = r/R 
weighted by the kinetic energy of the oscillation mode, 
and determines in which region the mode is trapped. 

- Column 10 - Z\ = l^J (x^) — {x)^, that takes values be- 
tween and 1, vanishes for a perfectly trapped mode. 

Since discontinuities in the sound speed derivatives in 
the stellar interior (such as those introduced by the bound- 
aries of the convective regions and by the second helium ion- 
ization zone) produce periodi c signatures in the frequencies 
of low degree modes (see e.g. lGoughllT990h . we decided to 
provide, as well, files containing the time variations of some 
quantities that may be of interest in the seismic analysis of 
stellar models. So, for each stellar parameter set (M,X,Z), 
but only for the MS-lifetime, we built a file whose name is 



ym- 



An individual file is provided for each stellar model whose 
name is formed with the stellar parameters, for example: ESTA/CoRoT" at |http://www.astro.up.pt/corot/ntools/info.html| 



Document "Description of the file formats used within 
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Fig. 6 Upper panel: equilibrium abundance of ^He for models with 
masses from 0.8 to 8 M^. Lower panel: equilibrium abundance of '''N 
(solid lines) and '^C (dotted lines) for each of the considered stellar 
masses. The thicker lines correspond to the chemical profile in a 2 Mq 
model. 



formed, as for HR diagram ones, with the values of the stel- 
lar parameters (for instance " ml.00Z0.02X0.70-ACC.txt"), 
and whose content in the columns 7-13 is the following: 

- Cutoff frequency (in ;uHz) at the photosphere (column 
7): Vac = c/{4nHp) where c is the sound speed, and Hp 
the pressure scale height. 

- Acoustic radius in seconds (column 8): r^c = j^\/cdr 
where R is the stellar radius at the photosphere. 



- Acoustic depth in seconds of the border of the convec- 
tive core (column 10) corresponding to the linear relative 
radius rcc/^ given in column 9. 

- Acoustic depth in seconds of the bottom of the convec- 
tive envelope (column 12) corresponding to the linear 
relative radius Tee//? given in column 11. 

- Second He ionization region: linear relative radius 
(rHeii/^) in column 12, and the corresponding acoustic 
depth ThcII in column 13. 

The four first lines, and the columns 1-6 are the same as in 
HRD2-files. 

In principle, only oscillation modes with frequencies lower 
than the acoustic cutoff one are trapped in the stellar interior 
and should be observed. As an indication of the order of 
magnitude of the expected oscillation frequencies we show 
in Fig. |7] the variation with time of the cutoff frequency for 
models with masses smaller than 2. Mq. 



8000 




0.6 0.4 0.2 



Fig. 7 Evolution of the cutoff frequency along the MS evolution of 
models with masses from 0.8 Mq (upper curve) to 2 (lower curve). 



5.1 Effect of atmospheric boundary conditions 

For the chemical composition X = 0.70 and Z = 0.02 there 
are two types of stellar models depending on whether the 
boundary conditions at T = Tgff are given by Kurucz at- 
mosphere models or by grey ones (with Eddington's law). 
We showed in Fig. [T] and |2] that the effect on the HRD lo- 
cation was significant only for T^q ^ 6300 K. The effect on 
the oscillation frequencies may be, however, important even 
for higher temperatures. As an example we show in Fig. [8] 
the difference between oscillation frequencies for 1.8 Mq 
models in mid-MS. The difference in the global parameters 
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for these Eddington and Kurucz stellar models are AR/R = 
4 X 10-"^, AL/L = 1.4 X 10"^ and AX^ < 10"^. Neverthe- 
less, the density atT — T^fi provided by the two atmosphere 
models is different enough to affect a significant fraction of 
the star, leading to an important variation of the oscillation 
frequencies. Note that for this comparison only the subpho- 
tospheric stellar structure was considered in the computa- 
tion of the oscillation frequencies. The differences at high 
frequency will be even larger if the atmosphere structure is 
included in the calculation. 




200 400 600 800 1000 1200 

i/{/xHz) 



Fig. 8 Difference between the oscillation frequencies of 1.8 Mq mod- 
els computed with Eddington's grey law (Ve) and with Kurucz's atmo- 
spheres (Vk). 



6 Conclusions 

In order to provide a reference grid of stellar models for the 
first interpretation of the CoRoT data, we have computed 
three sets of grey stellar models for 29 different masses from 
0.8 to 8 Mq. Each grid is computed with a different chemi- 
cal composition, and altogether they span metallicity values 
from [M/H]=-0.45 to -1-0.07 (taking as reference the value of 
Z/X=0.0245 for the Sun given by GN93). There are 3 x 29 x 
3 = 261 files containing data on time evolution of global 
quantities, around 3050 files with the detailed stellar struc- 
ture (interior and atmosphere), and other 3050 with the cor- 
responding oscillation frequencies for low degree modes (£ = 
0, 1,2, and 3). In addition, a set of stellar models with the 
Kurucz's atmospheres as boundary conditions were also com- 
puted, for the same values of stellar mass, but for only one 
chemical composition. That means additional ~2087 files, 
87 for the time evolution, and the rest of them for the stellar 
structure and the corresponding seismic properties. 
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